Transition to complete synchronization in phase coupled oscillators with nearest 
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We investigate synchronization in a Kuramoto-Iike model with nearest neighbour coupling. Upon 
analyzing the behaviour of individual oscillators at the onset of complete synchronization, we show 
that the time interval between bursts in the time dependence of the frequencies of the oscillators 
exhibits universal scaling and blows up at the critical coupling strength. We also bring out a key 
mechanism that leads to phase locking. Finally, we deduce forms for the phases and frequencies at 
the onset of complete synchronization. 
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Weakly coupled oscillators play an important 
role in understanding collective behaviour of large 
populations. They are often used to model the 
dynamics of a variety of systems that arise in na- 
ture, even though they are quite different. Syn- 
chronization is one of the interesting phenomena 
observed in these systems where the interacting 
oscillators under the influence of coupling would 
have a common frequency. Particularly, these 
systems show an extremely complex clustering 
behavior as a function of the coupling strength. 
In spite of their differences, the above mentioned 
systems can be described using simple models of 
coupled phase equations such as the Kuramoto 
model. This paper analyzes the behaviour of in- 
dividual oscillators in the vicinity of the critical 
coupling where all the oscillators evolve in syn- 
chrony with each other. 



I. INTRODUCTION 

Systems of coupled oscillators can describe problems in 
physics, chemistry, biology, neuroscience and other disci- 
plines. They have been widely used to model several phe- 
nomena such as: Josephson junction arrays, multimode 
lasers, vortex dynamics in fluids, biological information 
processes, neurodynamicsii^i^. These systems have been 
observed to synchronize themselves to a common fre- 
quency, when the coupling strength between these os- 
cillators is increased^'^'^'^'^ ' 1 1 1 2 1 1 3 ^ r^^ie synchroniza- 
tion features of many of the above mentioned systems, in 
spite of the diversity of the dynamics, might be described 
using simple models of weakly coupled phase oscillators 
such as the Kuramoto model^ii^. 

Finite range interactions are more realistic for the de- 



scription of many physical systems, although finite range 
coupled systems are difficult to analyze and to solve an- 
alytically. However, in order to figure out the collective 
phenomena when finite range interactions are considered, 
it is of fundamental importance to study and to under- 
stand the nearest neighbour interactions, which is the 
simplest form of the local interactions. In this context, 
a simplified version of the Kuramoto model with nearest 
neighbour coupling in a ring topology, which we shall re- 
fer to as locally coupled Kuramoto model (LCKM), repre- 
sents a good candidate to describe the dynamics of cou- 
pled systems with local interactions such as Josephson 
junctions, coupled lasers, neurons, chains with disorders, 
multi-cellular systems in biology and in communication 
systems^ii^ii^ii^. For instance, it has been shown that 
the equations of the resistively shunted junction which 
describe a ladder array of overdamped, critical-current 
disordered Josephson junctions that are current biased 
along the rungs of the ladder can be expressed by a 
LCKMii. In nearest neighbours coupled Rossler oscil- 
lators the phase synchronization can be described by the 
LCKMi^. Therefore, LCKM can provide a way to under- 
stand phase synchronization in coupled systems, for ex- 
ample, in locally coupled laser a^^'^° , where local interac- 
tions are dominant. Coupled phase oscillators described 
by LCKM can also be used to model the occurrence of 
travelling waves in neurons^i^. In communication sys- 
tems, unidirectionally coupled Kuramoto model can be 
used to describe an antenna array^^. Such unidirection- 
ally coupled Kuramoto models can be considered as a 
special case of the LCKM and it often mimics the same 
behaviour. 

One of the important features of the local model is 
that the properties of individual oscillators can be easily 
analyzed to study the collective dynamics while one has 
to rely on the average quantities, in a mean field approxi- 
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mation or by means of an order parameter, etc., as in the 
case of the usual Kuramoto model of long range interac- 
tions. Therefore, due to the difficulty in applying stan- 
dard techniques of statistical mechanics, one should look 
for a simple approach to understand the coupled system 
with local interactions by means of numerical study of 
a temporal behaviour of the individual oscillators. Such 
analysis is necessary in order to get a close picture on 
the effect of the local interactions on synchronization. In 
this case, numerical investigations can assist to figure out 
the mechanism of interactions at the stage of complete 
synchronization which in turns help to get an analytic 
solution. Earlier studies on the LCKM show several in- 
teresting features including tree structures with synchro- 
nized clusters, phase slips, bursting behaviour and sad- 
dle node bifurcation and so o n^'^i^^i^^ . There have been 
studies showing that neighbouring elements share dom- 
inating frequencies in their time spectra, and that this 
feature plays an important role in the dynamics of for- 
mation of clusters in the local model^^i^l. It has been 
found that the order parameter, which measures the evo- 
lution of the phases of the nearest neighbour oscillators, 
becomes maximum at the partial synchronization points 
inside the tree of synchronization^^. Very recently we 
developed a scheme based on the method of Lagrange 
multipliers to estimate the critical coupling strength for 
complete synchronization in the local Kuramoto model 
with different boundary conditions'^. 

In this paper we address the mechanism that leads to 
a complete synchronization in the Kuramoto model with 
local coupling. This is done by analyzing the behaviour 
of each individual oscillators at the onset of synchroniza- 
tion. For this purpose we consider the equations gov- 
erning the phase differences at the onset of synchroniza- 
tion. In particular, we identify that the cosine of only 
one among the phase differences becomes zero. Based on 
this property we derive the expression for the time in- 
terval between bursting behaviour of the instantaneous 
frequencies of each individual oscillators in the vicinity 
of critical coupling strength. Our analysis shows that 
the transition to complete synchronization occurs due to 
a saddle node bifurcation in agreement with the earlier 
studies. Further we deduce the expressions for the phases 
and frequencies of the individual oscillators at the onset 
of complete synchronization. 

This paper is organized as follows. In Sec. |ll] we 
present a brief overview on the dynamics of local Ku- 
ramoto model. Then wc analyze the behaviour of the 
phase differences and the time interval between succes- 
sive bursts at the transition to complete synchronization. 
In particular, we point out the mechanism that lead to 
complete phase locking at the critical coupling strength. 
Based on this we deduce the forms of phases and frequen- 
cies at the onset of synchronization. Finally, in Sec. IIIII 
we give a summary of the results and conclusions. 



II. BEHAVIOUR OF PHASES AND 
FREQUENCIES AT THE ONSET OF 
SYNCHRONIZATION 

Even when there has been an extensive exploration of 
the dynamics of the Kuramoto model (global coupling 
among all oscillators), the local model of nearest neigh- 
bour interactions, which can be considered as a diffu- 
sive version of the Kuramoto model, has been receiv- 
ing attention only recently. The LCKM is expressed 

j^ ^25.26.27.28.29 . 

6, =Lo, + y [sin(0,+i - e,) + sin(ft,_i - B,)], (1) 

here uji are the natural frequencies, K is the coupling 
strength, 9i is the instantaneous phase, 9i is the instan- 
taneous frequency and i = 1,2, . . . , N . Many interest- 
ing features of the LCKM remain unknown, especially 
an analytic solutionis, which would be of great impor- 
tance in understanding the mechanism that leads to syn- 
chronization. In order to find such an analytic solution, 
one should study carefully the temporal evolution of fre- 
quency and phase of each individual oscillator in the 
neighbourhood of the critical coupling for complete syn- 
chronization. 

If we consider the oscillators in a ring, with periodic 
boundary conditions ^i+Ar = 9i, the nonidentical oscil- 
lators H]) cluster in time averaged frequency, until they 
completely synchronize to a common value of the aver- 
age frequency ojn = jfJ2^=i'^i^ ^ critical coupling 
^^j25^^^^ At K > Kc the phases and the frequen- 



3 



V 

\ 


iV = 15 ■ 

Kc ■ 
(14-6)/ 




(7-13) 



K 

FIG. 1: (Color online) Synchronization tree for a system of 
15 oscillators. 

cics are time independent and all the oscillators remain 
synchronized. In Fig.[l] we show the synchronization tree 
for a periodic system with = 15 oscillators, where the 
elements which compose each one of the major clusters 
that merge into one at Kc, are indicated in each branch. 
In terms of phase differences (jji = 0^+1 — 6i, system 
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([T]), can be rewritten as: 
K 

(pi = cjj+i - cj, + y [sin(/)i„i 



2sin0i + sin(^i+i] , (2) 



with = at Kc for i = l,2,...,iV. In addi- 
tion, all quantities (j)*^ 0* and 6*, which become time 
independen t^^i^^i^^i^^ at the critical coupling, remain like 
that k > Kc when 4>* = and 9* = wq- Earlier attempts 
to get a solution of the above eq. ([2]) show that for only 
two oscillators which have phase difference 0^ = ^;*— (^jLi, 
results I sin 0^ | = 1— and indeed this is a necessary con- 
dition for eq. ^ to have a phase-locked solution. This 
fact has been used by Daniels et alii to estimate the 
value of critical coupling strength ifc, at which the tran- 
sition to complete synchronization occurs. However the 
determination of which two oscillators among TV oscil- 
lators that have |sin0^| ~ 1, remains difficult. From 
the study of the temporal evolution of phases and fre- 
quencies of each individual oscillator, it has also been 
found numerically that, at the onset of synchronization 
K < Kc, the values of 9i{t) and 4'i{t) remain equal to 
Wo and zero, respectively, for a certain time interval T. 
During this time T a stable phase-locked solution ex- 
ists, then they bursti^i^i^, and this stable phase-locked 
solution is lost. In between bursts, the phases remain 
fixed and then they have an abrupt change (phase-slip 
behaviour) by an amount which depends on the initial 
values of the frequencies ujr^-^, corresponding to the 
burst in the frequencies, while the quantities 0i = 

and 0i = are always preserved by the topol- 

ogy. Integrated with the above information, it has been 
shown by numerical investigation that the time interval 
T blows up as K becomes close to and T ^ oo at 
Kc- All these information leads to conclude that there is 
a saddle-node bifurcation at Kc and the synchronization- 
desynchronization transition at the critical coupling can 
be interpreted using this knowledge. 

In this work, we perform numerical investigations of 
the temporal evolution of the phases and frequencies for 
the individual oscillators in order to arrive to specific 
conditions which will lead to criteria to obtain an analytic 
solution. A detailed study of all quantities sine/)* at Kc 
for several values of TV and for different sets of Ui, shows 
that there is only one value of phase difference between 
two neighbouring oscillators (f)*i = 6*;*^^ — 9^ for which 
I sin (/)^ I = 1, while for all other values, i ^ I |sin(/)*| ^ 
1. In Fig. [21 we show sin^^ for a case of TV = 15 as 
time progresses at the critical coupling Kc, with the same 
initial frequencies of Fig. 1. We see that the value of 
I sin 0;* I = 1 , is for Z = 13 and that this quantity | sin 0j* | = 
1 holds for only one value of phase difference 0/ = 7r/2 
where these two oscillators I -\- 1 and I belong to different 
clusters, and these two nearest neighbours oscillators are 
always at the borders between the major clusters that 
merge at K^ which can be seen from Fig. 1. We find 
the same result for different initial frequencies oji and for 
different values of TV. In addition, the sign of sin^j* is 
negative for uji > uji^i and positive for the reverse. 



1 

0.5 


-0.5 
-1 



sin (pi = sin i 



>13 



A^= 15 



10 



20 30 
time xlO^ 



40 



50 



FIG. 2: (Color online) Values of sin (pi a.t K > Kc for a system 
of 15 oscillators [see Fig. [T] . 



The knowledge of the burst and phase slip (in the vicin- 
ity of Kc) of the quantities 4>i{t) and (piit), respectively, 
as well as the finding of | sin(/);*| = 1 (at Kc), will allow 
us to rewrite equation ([2|), for the index I as: 



B (A -2 sin </),), 



(3) 



where A 



3|(u;,+^-c^,)l +sin0;„i -fsin0;+i and B = f . 
Eq. ^ takes the form of a phase synchronization of two 
coupled limit-cycles^. At Kc, 0;* = ±7r/2 and, (j>i_i and 
(pi+i are constants and time independent and A = 2. A 
detailed numerical study shows also that, at the onset 
of synchronization, A 2 and the values of (pi, (pi-i 
and (pi-\-i remain equal to their values at Kc, for a time 
interval T. The values of Kc and A for different number 
of oscillators TV from numerical simulations are tabulated 
in Table, m It is clear that in all the cases, A ~2 when K 



TABLE I: Calculated values of Kc and A for different values 
of TV. 



N 


Kc 


A 


N 


Kc 


A 


3 


0.85041227 


1.9994 


20 


4.95830014 


2.0002 


5 


3.17082713 


2.0001 


25 


3.64106038 


1.9989 


10 


3.54701035 


1.9996 


50 


9.45720049 


1.9993 


15 


3.87023866 


2.0000 


100 


12.7232087 


1.9985 



approaches Kc- The relation A ~ 2 is found to be valid 
for different choices of initial frequencies LOi for each TV 
in the vicinity of Kc- Further, it should be noted that 
when the time interval T —f (X) one can find that A = 2. 
The time interval T can be found analytically, according 
to eq. (3), to be 
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r~ o/iv^ (4) 
kVaVa^ 

In Fig. [3] we clearly see that T blows up as A becomes 
close to 2 (where K goes to Kc), for the case of TV = 
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15. We find that T blows up as [A — 2) which is a 
numerical proof that a saddle-node bifurcation occurs at 
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FIG. 3: (Color online) logT versus log {A — 2), which shows 
the divergence of the time interval T when A approaches 2 
with a slope ~ —0.5 . 

Kc- Assuming that sin(/);_i and sin(/)i+i remain constant 
for a time interval T, in the vicinity of Kc, and equal to 
their values at Kc (which has been verified numerically) , 
we find that w Kc- Table HTl shows this fact where 



TABLE II: Calculated values of ^ for Af = 15 oscillators at 
the vicinity of = 3.870238658. 



K 


A 


AK 
2 


\K, - A'l 


3.869480136 


2.0007000 


3.870834454 


5.960x10" 


4 


3.870198414 


2.0000350 


3.870266142 


2.750x10" 


5 


3.870226709 


2.0000100 


3.870246060 


7.402x10" 


6 


3.870237122 


2.0000010 


3.870238909 


2.480x10" 


7 


3.870238325 


2.0000003 


3.870238727 


6.920x10" 
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the error is small and decreases as K approaches Kc 
Therefore, eq. @ takes the form 



^/2^c^Kc-K' 



(5) 



So that, within a good approximation, the periodic time 
interval T blows up as [Kc — K)~^-^, in good agreement 
with the numerical calculation by Zheng et. al^ii^, show- 
ing that a saddle-node bifurcation occurs at K^. 
Therefore, eq. ^ can be written as 

<Pi^ L{Kc- K sin (6) 

which can be solved analytically and its solution reads 

!tan (5 ait) ± K 



2 arctan 



Kc 



(7a) 



and 



L 



o? sec^ (iaLt) 



Kc 



[Q;tan(iQ!Lt) ii^jj' 



(7b) 



where a ~ \J K^ — K"^ and = § 
show that, at K^ the values sin (6 



Eqs. dZa]) and ([7b)) 
= ±1 which lead to 
0;* = 0. Also, it can be seen that in the vicinity of K^ 
sin(/!)/ = ±1 and 0/ = for a period T. The (-I-) sign in 
Eqs. (7) is corresponding to the case wj+i > and the 
(— ) sign for the reverse. 

In order to understand the mechanism of full syn- 
chronization which occurs at Kc, we use the fact that 
sint^j* = ±1 and each Q* ~ i^Q, where these quantities re- 
main unchanged for T in the vicinity of Kc- Hence, from 
system (1), we are able to get the following relations: 



N-l 



Kc 



(8a) 



m—l 
„ l-l 



— ^(cjo - w;_„_i) ± sin0*. (8b) 



Using this fact, we write the following equations, in ad- 
dition to eq. ([7]): 



(Pl-n 
4>l-n 
h+m 



sin ^(a„±sin0i). 



cos (pi (pi 



yr^(a7±siii^I)^' 
sin"-^(a„i ± sin (pi), 

cos (pi (pi 
^1 - [am ± sin(/)/)2 



(9a) 
(9b) 
(9c) 
(9d) 



where a„ 
1 and am 



-3 



^^(wo-w/_j_i) with n = 1, 2,3, . . . , Z- 

7f;I]j1i(^o-w/+j)withm= 1,2,3, ... ,7V- 
l. It is clearly seen that according to the above equation, 
each (pi can be expressed in terms of (pi and consequently 
each (pi can be expressed in terms of (pi and (pi . There- 
fore, all values of (pi will be shifted from each other by 
some constant which is determined by the location of the 
indexes / — n and I + m relative to oscillators with in- 
dexes I and / + 1. This is shown in Fig. (2), where sm. (pi 
values are shifted from each others at Kc- Therefore, at 
Kc, what occurs to (pi and 0/ due to saddle-node bifur- 
cation diffuses through the ring via interaction between 
neighbouring oscillators. This means that, at the vicinity 
of Kc, the value of (pi has an abrupt change after being 
constant for a time T, caused by a burst behaviour of (pi 
after being zero for the same time interval T. The abrupt 
change in (pi produces a sudden change in the values of (pi 
of their neighbours, while the bursting behaviour of (pi in 
turn yields bursts in (i ^ I). In order to demonstrate 
this fact, we plot the temporal evolution of (pis and (pe, in 
the vicinity of Kc, according to numerical simulation of 
eq. (1) in Fig. 4a while we plot both quantities according 
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FIG. 4: (Color online) Time evolution of ^13 and i^e according 
to (a) system (1) and (b) equations (7b) and (9b), at -ft' = 
3.870226709, for 15 oscillators with the same initial conditions 
of Fig. 1 



to eq. (7b) and eq. (9b) in Fig. 4b. As shown in Fig. 4, 
the results of numerical simulation agrees with that from 
the analytic solution. The above mentioned behaviour 
is reflected in the time dependence of the 6[s, which in 
turn remain equal to loq for a time T and burst around 
uq corresponding to the burst of 0;. Henceforward, we 
argue that it is the behaviour of and 4>i, which drives 
the system to fall into full synchronization. 



III. SUMMARY AND CONCLUSIONS 

In summary, we have analyzed the conditions on the 
phase differences for the onset of complete synchroniza- 
tion at the critical coupling strength in a Kuramoto-like 
model with nearest neighbour coupling. Such condition, 
which is I sin (f>'^\ = 1 (or cos 0;* = 0) , allows us to solve 
the equations of the phase differences ^ analytically. 
Also, we found that full synchronization occurs always 
when the quantity A = 2 at Kc- Due to the diffusive 
nature of the LCKM, complete synchronization of all os- 



cillators to a common value can be interpreted and un- 
derstood once we have an analytic forms for and </);. 
However, it is still difficult to determine analytically the 
number of oscillators in each cluster which merge into 
one at Kc- Therefore, one cannot allocate straightfor- 
wardly the two nearest neighbour oscillators which would 
have I sin0;*| = 1. On the other hand, a detailed numer- 
ical study on the temporal evolution of phases, phase- 
differences and frequencies of oscillators at the borders of 
the clusters that merge into larger one at onset of com- 
plete synchronization helps us to determine the neigh- 
bouring oscillators which have sint/);* = ±1. Such anal- 
ysis can also be used to understand the partial synchro- 
nization that leads to the formation of small clusters for 
coupling strengths below the critical coupling strength 
Kc- Of course analysis of the simplest case of locally 
coupled phase oscillators can help to understand mod- 
els with local interactions where amplitudes and phases 
are include d^^i^^'^^i^° . In such cases, a detailed study 
of the time evolution of amplitudes and phases can re- 
veal a better understanding of the mechanism of synchro- 
nization. The present analysis can also be applicable 
to models in higher dimensions such as that for dislo- 
cations in solids which includes local nearest neighbour 
interactions^. Furthermore, the present approach can 
be extended to understand the underlying mechanism in 
the case of locally coupled Kuramoto models with time 
delay^ (or phase delay) introduced between the coupled 
oscillators. In addition, the mechanism of synchroniza- 
tion in LCKM for open and fixed boundaries can be stud- 
ied in a similar manner to the present work as well as for 
the case of unidirectional LCKM. We also want to men- 
tion that the scaling law given by eq. (5) has been found 
experimentally in a transition to phase synchronization 
in CO2 lasers^ and in electronic circuits^^i^i. On the 
other hand one can not make a direct comparison be- 
tween the mechanism of synchronization discussed here 
in LCKM and the scaling law that has been found in ex- 
periments since the physical systems are not necessarily 
the same. 
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